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Abstract 

We describe a way to optimize the chiral behavior of Wilson-type lat- 
tice fermion actions by studying the low energy real eigenmodes of the 
Dirac operator. We find a candidate action, the clover action with fat 
links with a tuned clover term. We present a calculation of spectroscopy 
and matrix elements at Wilson gauge coupling f3 = 5.7. The action shows 
good scaling behavior. 



1 Introduction 



There is an increasing accumulation of evidence from lattice simulations of the 
importance of the topological properties of the QCD vacuum @. The (lattice) 
pure gauge QCD vacuum is filled with instantons of average radius ~ 0.3 fm, 
with a density of about 1 fm , and they are responsible for a large part of 
chiral symmetry breaking in QCD ||, The extent to which the complete 
continuum phenomenology of instanton effects in hadronic physics is actually 
realized by QCD, and can be studied using lattice simulations, remains an open 
question, but its outline is present. 

Lattice artifacts associated with instantons also seem to be connected to 
some of the difficulties of numerical simulations with quenched QCD, namely, 
eigenmodes of the Dirac operator which occur away from zero bare quark mass, 
and which spoil the calculation of the fermion propagator, the so-called "excep- 
tional configurations" [Q, 

In this work we wish to look at topology and fermions on the lattice from a 
different perspective: given that instantons are responsible for chiral symmetry 
breaking, is it possible to optimize the lattice discretization of a fermion action 
with respect to its topological properties? Our goal is to shrink the range of 
bare quark mass over which the low lying real eigenmodes of the Dirac operator 
occur, on background gauge field configurations which are typical equilibrium 
configurations at some gauge coupling. 

The action we propose is the standard clover action, except that the gauge 
connections are replaced by APE-blocked || links, and the clover coefficient is 
tuned to optimize chiral properties. Specifically, 

S = y~^(m + 4)^(n)^(n) 

n 

- ~ 52$(n)(l - ^)V^ri)ip(n + y) + $(n)(l + l^V^n - ^{n - fi) 

C , 

+ tjj{n)i<T^ u P^(n)^(n) (1) 

with 

v<r\x)= {i-c)v!i n - x \x) 

+c/GY,^ fl (vS n - 1) (x)vl n - 1 \x + vW^ix + A)* 
+vS n - 1) (x - i>)^ n - 1] {x - i>)v^- 1] {x -0 + fi), (2) 
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with V^ n \x) projected back onto SU(3) and V^\n) = U^(n) the original link 
variable. We take c = 0.45 and N — 10 smearing steps, chosen because of our 
previous work in instantons ^]. This choice of parameters is not unique and 
might not even be optimal. P M „ is the usual clover set of links, but built of the 
V^'s. At Wilson gauge coupling 5.7 — 5.8, the best choice is C — 1.2, and it 
decreases to the tree-level C = 1 value at larger (3. We will also quote simulation 
data in terms of the hopping parameter, ^ = m + 4. 

The optimized action has a spread of low lying real eigenmodes with respect 
to the bare quark mass which is less then a third of the usual Wilson action. In 
terms of the square if the pion mass the spread is about three times smaller for 
the optimized action than for the standard Wilson action. The action has other 
good features, as well: the renormalization factors connecting lattice quantities 
to their continuum values appear to be very close to unity. The action also 
appears to require only about half the number of sparse matrix inversion steps 
as the usual clover action (at equivalent values of physical parameters). 

It is possible to tune the standard clover action (with the original links 
used as connections) to improve the interaction of quarks with instantons. This 
choice is not the best one. It happens that the best value of C for f3 < 6.0 
is large enough that a new class of exceptional configurations appear. These 
configurations compromise simulations. The connection between exceptional 
configurations and instantons is explored in detail in a companion work by us 
I- 

The action we propose is completely unimproved in its kinetic properties, 
so its dispersion relation and heavy quark mass artifacts are identical to those 
of the Wilson action. It would be very easy to improve it by beginning with a 
more complicated free fermion discretization. 

Fermion actions with fat links have received considerable attention in the 
past year. The first use of them we know of (although with a different mo- 
tivation) was in the simulations of QCD on cooled gauge fields by the MIT 
group fj. The MILC collaboration g, Or ginos and Toussaint || and Lagae 



and D. K. Sinclair) 10 have shown that modest fattening considerably improves 
flavor symmetry restoration in simulations with staggered fermions. We have 
used calculations of staggered spectroscopy on highly smoothed gauge config- 
urations to compare chiral symmetry breaking in SU (2) gauge theory and in 
instanton backgrounds [fxifl . Finally, all fixed point actions and approximate 
fixed point actions we know of [to , [l3]] for fermions seem to incorporate fat 
links. Fixed point fermions realize the index theorem and retain chiral symme- 
try at nonzero lattice spacing |5 , and so one way of viewing a fat link action 
is as an approximate FP action, which includes its chiral properties but not its 
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kinetic ones. 



An apparent drawback of fat link actions is the lack of a transfer matrix 
between consecutive time slices. N APE steps can mix links up to ±N timeslices 
away and a strict transfer matrix cannot be denned on time slices closer than 
2N lattice spacings (20 in our case). In practice we are concerned only with 
the exponential decay of correlation functions and we found asymptotic decay 
after 3-5 time slices in our simulation. This is expected if we realize that APE 
smearing is basically a random walk whose range can be estimated as \fNc w 
1.4. 

The outline of the paper is as follows: In Section 2 we review the continuum 
index theorem and describe how lattice fermions fail to reproduce it. We then 
describe how we find the real eigenmodes of the Wilson-Dirac operator. All 
this is basically a review, and experts may skip it. In Section 3 we describe 
the tests we performed to tune the action. Section 4 describes a calculation of 
spectroscopy and matrix elements at aT c = 1/4 (Wilson gauge coupling /3 = 5.7) 
using the new action. 



Measuring the Chiral Properties of Wilson- 
like Fermions 



The reader might recall that the local topological density 

Q( x ) — TTo o^fpc-F/w-Fper (3) 
OZIT 

is related to the divergence of the flavor singlet (with n / flavors) axial- vector 
current 

d^ij^ip = 2m-0i75-0 + 2in f q. (4) 
Integrating over all x, this relation implies a connection between the topological 
charge Q — J d i xq(x) and a mode sum, 

Q = mTr^ = ™Tr 75 1 = m V .{' 7 * /a (5) 
7 • D + m i\ s + m 

where f s are the eigenfunctions of the (antihermetian) Dirac operator 7 • D, 
7 • Df s = iX s f s , with A s real. The property {75,7 • D} = leads to the 
condition that fjjsfs = if 7^ 0, and if A s = we may choose f s to have a 
definite chirality, j 5 f s = ±f s . Thus it follows that 

Q= fhsf* =n+-n- (6) 

s,A s =0 
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where n+ and n_ are the number of zero eigenmodes with positive and negative 
chirality. This is the index theorem. 

On the lattice, essentially every statement in the preceding paragraph is 
contaminated by lattice artifacts fL5[ ], Here we focus on Wilson-like fermions, 
where chiral symmetry is broken by the addition of terms proportional to the 
Dirac scalar and/or tensor operators. The Wilson or clover fermion action 
analog of 7 • D, D W1 is neither Hermetian nor antihermetian and its eigenvalues 
are generally complex. D w can also have real eigenvalues. These real eigenvalues 
usually do not occur at zero bare quark mass. Their locations spread across a 
range of quark mass values. 

On smooth, isolated instanton background configurations, the location of the 
real eigenmode varies with the size p of the instanton. For large instantons, the 
low lying real mode occurs close to zero quark mass. Accompanying these near- 
zero modes are a set of modes which do not cluster around mo = 0, but around 
—arriQ ~ 0(1). These are "doubler modes." As one decreases the instanton size, 
the eigenmodes at arriQ ~ shift towards negative quark masses, approaching 
the doubler modes. This p-dependent mass shift and the location of the doubler 
modes all depend on the particular choice of lattice fermion action p|, [l6| |l7||. 
The shift of the low lying real eigenmode for a given instanton size for the Wilson 
action is larger than for the clover action, which has better chiral properties. 
Nevertheless, for both actions, as the instanton size decreases, sooner or later 
the low energy eigenmode shifts to large negative quark mass, approaches the 
doublers and eventually annihilates with one of the doublers - the fermion does 
not see the instanton any longer. 

On equilibrium background configurations, the shift of the eigenmode with p 
is accompanied by an overall /3-dependent mass shift. On top of that there is the 
spread of the eigenmodes according to the sizes of the background instantons. 
The spread of these low energy modes, if distinguishable from the doublers, 
characterize the amount of explicit chiral symmetry breaking of the fermionic 
action. On configurations with small lattice spacing the typical instanton is 
large in lattice units. The spread of the low energy modes is small and they 
are well separated from the doublers. The explicit chiral symmetry breaking 
of the fermionic action is small and controlled. On configurations with large 
lattice spacing the instantons are small (if present at all). The spread of the 
low lying modes is large, close to or overlapping with the doubler modes. The 
hadron spectrum on these configurations could numerically be similar to the 
continuum spectrum, but the physical mechanism behind it is very different 
from continuum QCD. 

Determination of the locations of the real roots of D w + m in a particular 
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background gauge configuration can be done very simply. We approximate 



P(m ) = (V'75'0> = Tr(fl„ + m ) *7s 



(7) 



with a noisy estimator: cast a random vector rji on each site i, compute (D w + 
m o)x = V an d measure 



as a function of m a . The eigenmode of D w is located at —m and its appear- 
ance is signaled by the appearance of a pole in P(mo). This is a variant on 
the standard method of measuring (tpip) in a lattice simulation. It has also 
seen considerable use by the Fermilab group Q in their studies of exceptional 
configurations (they use a flat source rji = constant, not a noisy source). The 
eigenfunction itself of a particular mode can be found by performing the inver- 
sion for a test mass very close to an eigenvalue. Then 



is saturated by the mode closest to the pole, and x'x is (proportional to) the 
probability density for that mode. 



3 Testing and Tuning Actions 

We study the locations of low energy real eigenmodes in equilibrium (quenched) 
gauge configurations. Here there are two issues: First, do the low lying modes 
separate from the doubler modes, and second, what is the spread of the low 
lying modes. 

The first question is important because, if the low lying modes do not sep- 
arate from the doubler modes, then the physics of chiral symmetry breaking 
with a lattice cutoff is different than in the continuum. The lattice theory is 
no longer just a crude approximation to the continuum, and one cannot speak 
of chiral symmetry breaking as being induced by instantons. Chiral symmetry 
is certainly broken in QCD for any value of the cutoff, including the strong 
coupling limit but the mechanism does not involve instantons. 

Note that if the low lying modes do not separate from the doubler modes, it 
does not make sense to talk about the spread of the low lying modes. 

In this paper we restrict ourselves to clover-like actions. Our goal is to tune 
the clover coefficient for good chiral behavior, i.e. to minimize the spread of the 




(8) 




(9) 
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low lying modes. We could attempt this program on the original configurations. 
However, it is known that rough gauge configurations with large clover coefficient 
are plagued by exceptional configurations. This problem is greatly reduced 
if the links of the configurations are smoothed by a series of APE smearing 
steps j5j. We also know from previous work jll| that an action with fat links 
is insensitive to short distance fluctuations, but still knows about instantons 
and the additional long distance behavior of the gauge field responsible for 
confinement. So we will begin with a fat link action of Eqn. [|. Throughout this 
paper we create the fat link by 10 APE smearing steps with smearing coefficient 
c = 0.45. 

Figure |l| shows the spread of the real eigenmodes at p = 5.8 as the function 
of the bare quark mass for clover coefficients C = 0.0 (Wilson action), C = 1.0, 
C = 1.2 and C = 1.4. The eigenmodes were located using the pseudoscalar 
density function of Eqn. || on 40 8 4 configurations. The horizontal ranges of the 
histograms are equal to the range the eigenmode search had been performed. 
To the right of the distributions is the confining phase. If one measures the 
pion mass in the positive mass region and extrapolates to zero with the 
bare mass mo, one finds that the bare mass at which m n = lies at a value 
m c located within the range of the low- lying real eigenmodes (j| . If a particular 
configuration has an eigenmode with mo > m c , and if one attempts to compute 
the quark propagator at m = mo, the propagator will be singular. This is 
(one kind of) exceptional configuration, and at a minimum it compromises the 
statistical averaging process inherent in the Monte Carlo simulation process. To 
the left of the distributions, in the negative quark mass region, are the doublers. 
In Figure |l| there is no indication of them, since the doublers did not show up 
within the range we scanned for eigenmodes. It appears that the doublers and 
low lying modes are well separated at P = 5.8. This seems to be in contradiction 
with the statement of Ref. |jl^| , where the authors claim that the real modes 
cover the whole investigated negative quark mass range. Of course, as there is 
a non-zero probability of finding an eigenmode anywhere, our claim is only that 
the distribution is strongly peaked at the low lying modes and the doublers are 
well separated. Since in Ref. |l9| no distribution histogram is published, we 
cannot tell if there is a real discrepancy between the results or it is the question 
of interpreting the data. One should also note that in |lj| thin link actions were 
used, so we are really discussing two different actions. We will return to this 
question in ||. 

To quantify the spread of the low energy real modes, in figure ||a we plot the 
average real eigenmode location as the function of the clover coefficient. The 
error bars here are not errors, they are the spread of the modes in m. Figure 
shows only the spread as the function of the clover coefficient. m c , where the 
pion becomes massless, lies somewhere in the middle of the mode distribution, 
closer to its right (large mass) end. The average eigenmode locations and the 
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upper end of the error bars in figure ga bracket m c . Smearing the link removes 
most of the additive mass renormalization even for the Wilson action, reducing 
m c ~ —0.95 for the thin link action to m c ~ —0.22 for our case. Adding 
a clover term to the action further reduces the additive mass renormalization, 
and even larger clover terms induce a positive mass renormalization. At j3 = 5.8 
the additive mass renormalization is minimal for c ss 1.2. This is also the value 
where the spread of the eigenmodes is minimal. To conclude, the results indicate 
that at (3 = 5.8 (lattice spacing a ~ 0.15 fm) the low lying and doubler modes are 
well separated and explicit chiral symmetry breaking is minimized with clover 
coefficient C — 1.2 with our fat link action. 

The situation is less convincing at (3 — 5.7 (lattice spacing a ~ 0.20 fm). 
Figure ^, again based on 40 8 4 configurations, shows the distribution of the 
eigenmodes at (3 = 5.7, clover coefficient C = 1.0, 1.2 and 1.4. Here the low lying 
modes are not as well separated than for (3 = 5.8. For C — 1.0 the distribution 
has a large tail extending towards negative quark masses. The situation is a 
bit better for C = 1.2. Since the low lying modes are not well separated from 
the doublers, it makes no sense to calculate the spread of the distribution. The 
continuum description of chiral symmetry breaking is about to break down at 
(3 = 5.7. If any of the distributions of figure || describes continuum physics, it 
is C = 1.2 where the overall mass renormalization is close to zero and the low 
lying modes are best separated form the doublers. 

Finally, at j3 = 5.55 (lattice spacing a ~ 0.24 fm) Fig. |^ shows that the 
distribution of eigenmodes is broad and the low lying modes and doublers are 
completely mixed. 

These figures show that it is not possible to make the lattice spacing greater 
than about 0.2 fm, and still retain the continuum-like description of chiral sym- 
metry breaking using a clover-like action. Other, better tuned actions might 
perform better at large lattice spacing. Our exploratory studies with a hyper- 
cubic fixed point jl3| action showed that the FP action is not better in terms 
of chiral symmetry breaking than the fat link clover action with C = 1.0. Its 
dispersion relation was improved at lattice spacing a — 0.36 fm, and it (and 
the clover action) showed only a small amount of scale violation in hyperfine 
splittings at large lattice spacing. This just shows that scaling or near scaling of 
a few quantities does not guarantee that the physical mechanism responsible for 
chiral symmetry breaking is the same as in the continuum. This should serve 
as a warning sign for any calculations at large lattice spacing with actions of 
untested chiral properties. 

In principle, one should optimize not the spread of the bare quark mass, but 
the spread of eigenmodes in terms of some physical observable, like the pion 
mass. However, at large fattening, the relation between the bare quark mass 
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and the pion mass shows little variation with C. 

The usual nonperturbative tuning of the clover term p0| attempts to opti- 
mize current algebra via PCAC relations. It is not clear to us how this is related 
to what we do, since it has been done so far only for thin link actions. We know 
that the optimal C by our criterion for a thin link action is larger than the 
usual nonperturbative tuning, for example 2.5 vs 2.0 at /3 = 5.8. Our criterion 
is based on a direct attack on low energy real eigenmodes. To the extent that 
the dynamics of QCD for light quark masses is dominated by this physics, we 
feel our optimization criterion is well founded. 

We do not know if there is a connection between this tuning method and 
the Ginsparg- Wilson |2l| realization of chiral symmetry on the lattice. The free 
field limit of these fat link actions is identical to the standard Wilson action, 
and so we are doing something different from improving the chiral properties of 
the interacting theory by improving the chiral properties of the free theory. 



4 Scaling Tests 

We have done a calculation of spectroscopy and matrix elements at lattice spac- 
ing aT c = 1/4 [fi = 5.7 with the Wilson gauge action) using the fat link clover 
action with C = 1.2. This is on the edge of the allowed range of lattice spacings 
according to our criterion of the last section. By itself, one lattice spacing is not 
a scaling test, but we can combine our results with those from other actions to 
compare the new action to them. 



4.1 Spectroscopy 

The spectroscopy measurement is entirely straightforward. The lattice vol- 
ume was 8 3 x 24. The data set was 80 lattices. We gauge fixed to Coulomb 
gauge and used a Gaussian independent particle source wave function ip(r) = 
exp(— (r/ro) 2 ) with ro = 2. We used pointlike sinks projected onto low mo- 
mentum states. We used naive currents (^75^, etc.) for interpolating fields. 
The spectra appeared to be asymptotic (as shown by good (correlated) fits to a 
single exponential) beginning at t ~ 3 — 5 and the best fits were selected using 
the HEMCGC criterion ||. 

Our fiducials for comparison are Wilson action and clover action quenched 
spectroscopy. We have tried to restrict the data we used for comparison to 
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Figure 1: Locations of the real eigenmodes at /? = 5.8 of fat link Wilson action 
(a); fat link clover fermions C=1.0 (b); C = 1.2 (c); C = 1.4 (d). 
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lattices with the proper physical volume. 

We can roughly estimate the critical bare quark mass (at which the pion is 
massless) by linearly extrapolating to zero in mo- We also estimated the 
critical bare mass using the PCAC relation 

V p • (^7 5 -0(O)V'757 '^(x)) = 2m q (ip-f 5 'ip(0)'iJ)^ 5 tp(x)) . (10) 

which, going to the lattice and following |23) is done by fitting the pseudoscalar 
source-pseudoscalar sink to 

P(t) = Z(exp(-TO w t) + exp(-m n (N t - t))) (11) 

and the pseudoscalar source-axial sink to 

= |P^£z(exp(-(TM) - exp(-m T (iV t - t))) (12) 

to extract m q . Fig. ^| shows the squared pion mass vs bare quark mass. The 
quark masses from the local pseudoscalar and axial currents are also shown. 

Fits of the squared pion mass (am,) 2 = B(l/n — 1/k c ) give B = 1.55(2) 
and k c = 0.12515(7). The quark mass is fit to am q — A(1/k — 1/k c ) an d 
we find A — 0.447(7), k c = 0.12515(6). In free field theory, we would expect 
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Figure 3: Locations of the real eigenmodes at j3 — 5.7 of fat link clover fermions 
C=1.0 (a); C= 1.2 (b); C = 1.4 (c). 
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Figure 4: Locations of the real cigenmodes at [3 = 5.55 of fat link clover fcrmions 
C=1.0 (a); C= 1.2 (b). 
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Figure 5: Bare squared pion mass (crosses) and quark mass from Eq. 12 (dia- 
monds) vs hopping parameter for the C=1.2 action. The octagon shows n c . 
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CD * 
CD □ 





1 3 3 4 5 

m n / T c 

Figure 6: Octagons show m p /T c vs. m w /T c for the C=1.2 action at aT c = 1/4. 
Also shown are Wilson action data, with crosses for aT c = 1/4, squares for 
aT c = 1/8, and fancy crosses for aT c = 1/12. 



A = 1/2, k c = 1/8, and this is our first hint that perturbative corrections to 
bare quantities are small. 

As a scaling test we compare m p /T c and /T c vs. m^/Tc for the C=1.2 
action in Figs. | and H 

We compare scaling violations in hyperfine splittings by interpolating our 
data to fixed it/ p mass ratios and plotting the N/p mass ratio vs. m p a. We do 
this at three tt/ p mass ratios, 0.80, 0.70 and 0.60, in Fig. ||. In these figures the 
diamonds are Wilson action data in lattices of fixed physical size, 8 3 at (3 = 5.7 
16 3 at /3 = 6.0 [|D 24 3 at ft = 6.3 p6| , and the crosses are data in various 
larger lattices: 16 3 and 24 3 at = 5.7 and 32 3 at j3 = 6.17 Q, 24 3 at (3 = 6.0 
[ p5[ . When they are present the data points from larger lattices illustrate the 
danger of performing scaling tests with data from different volumes. The bursts 
are from the nonperturbatively improved clover action of Refs. J27j] and J28|. 
The squares show the C=1.2 action. 

We give a table of masses from the C=1.2 action in Table [l]. 

We conclude that the fat link clover action has at least as good scaling 
behavior as the usual nonperturbatively tuned clover action. 
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iJJ -7T/ 1 C 

Figure 7: tun /T c vs. m^/T c for the C=1.2 action and the Wilson action, 
labelled as in Fig. 6. 



We also noticed that the new action seems to be more convergent than the 
thin link clover action: the same biconjugate gradient code needs about half as 
many steps to converge to the same residue as the usual thin link clover action, 
for the same ir/p mass ratio. 

We did not encounter any exceptional configurations with the fat link action 
for m^/nip > 0.56. 



4.2 Dispersion Relations 

To view the dispersion relation, we first plot E(p), the energy of the state 
produced with spatial momentum p, as a function of \p\. The result for the 
C=1.2 action at k — 0.118 is compared to the free dispersion relation at aT c = 
1/4 in Fig. ||. It looks very similar to results from the Wilson action at the 
same lattice spacing-an entirely unsurprising result. 

This behavior is quantified by measuring the squared speed of light, c 2 = 
(E(p) 2 —m 2 ) /p 2 , for p = (1, 0, 0). We do this by performing a correlated fit to the 
two propagators. The result is presented in Fig. |l^ and shows the worsening of 
the dispersion relation at larger quark mass, characteristic Wilson/clover kinetic 
behavior. 
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Figure 8: A scaling test for the C=1.2 action (square) vs. Wilson actions on 
lattices of fixed physical size (diamonds) and larger volumes (crosses), and the 
nonperturbatively improved clover action (bursts). Data are interpolated to 
n/p = 0.80 (a), 0.70 (b), and 0.60 (c). 
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Figure 9: Dispersion relation for hadrons at aT c = 1/4 (a ~ 0.18 fm) from 
the C—1.2 action. The curves are the continuum dispersion relation for the 
appropriate (measured) hadron mass. 



K 


PS 


V 


N 


A 


0.118 
0.120 
0.122 
0.123 


0.868( 3) 
0.730( 4) 
0.567( 5) 
0.467( 6) 


1.053( 6) 
0.965( 8) 
0.880(12) 
0.837(17) 


1.589(13) 
1.415(19) 
1.219(28) 
1.093(36) 


1.725(19) 
1.593(22) 
1.501(24) 
1.432(27) 



Table 1: Table of best-fit masses, C = 1.2, f3 = 5.7. 
4.3 Renormalization factors 

We now turn to a (rather naive) set of measurements of simple matrix elements. 

Formally, the fat link action can be regarded as "just another" 0(a 2 ) im- 
proved action, if one assumes that the coefficient of the clover term is actually 
C = 1 + 0(g 2 ). In perturbation theory, which involves the vector potential 

rather than the link, the action shows the usual cancellation between the 
clover term and the scalar -ipipA term. The spectrum is 0(a 2 ); matrix elements 
using "rotated fields" ip — > (1 — (7 • D)/2)ip are 0(a 2 ) improved. Of course, the 
lattice-to-continuum renormalization Z-factors are different than the usual thin 
link clover Z-factors. 

In principle, one could calculate all Z-factors using lattice perturbation the- 
ory. The Feynman rules differ from the rules for the usual clover action, only 
in that the vertices are multiplied by form factors, which are functions of the 
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Figure 10: Squared speed of light vs. hadron mass in units of T c , for (a) 
pseudoscalars, (b) vectors) and (c) protons, from the C=1.2 action at /? = 5.7. 
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gluonic momentum. We have not tried to do this yet. Since we are working at 
fairly strong coupling, we felt that it would be easier to compute the Z-factors 
nonpcrturbatively, beginning with the vector current and then computing the 
axial renormalization factor using Ward identities. Since this is just the first 
calculation using this action, we will restrict ourselves to the naive (local) cur- 
rents. 



We begin with the vector current. The conserved (Noether) current for the 
fat link action is the same as the Wilson current, just with fat links instead of 
thin links, 

J c r S ( n ) = ^W(n)(7M-l)^W(« + A)-^ + «(7, + l)^WW). (13) 

This current is conserved but not improved. We choose to define the vector 
current Z-factor from the ratio of forward matrix elements 

V (n(T)J (t)n(0)) { > 

for < t < T. The local current is defined as J l u°(n) — ij)(ri)'Y l j l ip(n). We 
measured the Z-factors using 20 configurations on a periodic 8 3 x 24 lattice 
with T = 10 and averaged 9 > t > 1. The results of a jacknife analysis are 
shown in Table ^. 
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The fact that Z cons ^ 1 is a finite size effect. In the numerator of Eqn 
both quarks can propagate in both time directions around the torus, but the 
forwards and backwards paths contribute differently to the denominator. We 
correct this phenomenologically by computing the ratio z loc /Z cons and display 
it in the table. 



Lattice perturbation theory which correctly takes into account the residue 
of the pole of the massive quark [^9| predicts that this ratio is equal to 

Zv = ^-^Z V (15) 

where the (1 — 6k) factor is the residue (recall k c = .125) and Zy = l + ag 2 + . . .. 
As Table g shows, the correction to the tree level formula, Z, which is what is 
usually quoted as the Z-factor, differs from unity by less than than two per cent 
at k = 0.123. 



We also attempted to measure the axial vector current renormalization using 
the chiral Ward identity |3(J. With 50 lattices, we only had a signal at k = 0.118, 
where we found Za = 1.210(8). This is again quite close to the pure kinetic 
result of 1.237. Dividing out the kinematic factors, this would give Za = 0.978. 
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K 


^cons 


Z loc 




i (j« 

2 k 


0.118 


1.020(1) 


1.253(2) 


1.228 


1.237 


0.120 


1.032(2) 


1.185(4) 


1.148 


1.166 


0.122 


1.054(5) 


1.129(6) 


1.071 


1.098 


0.123 


1.071(6) 


1.127(8) 


1.052 


1.065 



Table 2: Vector current renormalization factors, C = 1.2, (3 = 5.7. 

At j3 = 6.0, Ref. |l| found Z v = 0.824, Z A = 1.09 for the usual clover 
action, and improved operators. 

It would be very interesting to measure the mixing of left- and right-handed 
operators (for Br-) in this action. 



5 Conclusions 

We have shown that via a combination of fattening the links and tuning the 
magnitude of the clover term, it is possible to optimize the chiral behavior of 
the clover lattice fermion action. The example we presented reduces the spread 
in the low energy real eigenmodes by about a factor of three in units of the 
squared pion mass, compared to the usual Wilson action. We fixed the fattening 
and varied the clover coefficient, but it is clear that a real optimization would 
involve varying both factors. 

The same procedure can be applied to the usual clover action, without fat 
links, but the value of the clover term C which minimizes the spread of real 
eigenmodes due to instantons is so large that new kinds of configurations, whose 
singular behavior is unrelated to fermion modes sitting on individual instantons, 
become exceptional ||. Thus this choice of action does not represent a good 
choice for improvement. 

The method of construction exposes an apparent fundamental upper limit 
for the lattice spacing in a QCD simulation of about 0.2 fm. At larger lattice 
spacing, the doubler modes and the low energy modes do not show a clean 
separation, and the mechanism of chiral symmetry breaking is qualitatively 
different than in the continuum. The Wilson and clover actions, and the one 
improved action we tested, the hypercubic action of Ref. (lj| , failed at lattice 
spacing 0.24 fm. This test does not apply to actions involving heavy quarks, 
for which chiral symmetry is (presumably) not important, but before any other 
improved action can be said to reproduce continuum light hadron physics at 
some lattice spacing, it should show a separation between the spectrum of its 
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near low energy modes and its doubler modes. 

The "fat clover" action we have presented in this paper, as an example of 
an action with improved chiral properties, has many other nice properties as 
well. It appears to exhibit scaling of hyperfine splittings at 0.2 fm. The action 
is quite insensitive to the UV behavior of the underlying gauge field. It has 
Z-factors for simple matrix elements which are very close to unity. And finally, 
the amount of resources required to construct propagators is halved compared 
to the usual clover action, at equivalent parameter values. It would be trivial 
to modify any existing clover code to improve it in the way we have described. 
Of course, its kinetic properties, including the dispersion relation, artifacts at 
large am q , and power law scaling violations in matrix elements, are unchanged 
from the standard clover action's. 
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